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Abstract. In this paper we explain that the natural filtration of a continuous 
Hunt process is continuous, and show that martingales over such a filtration are 
continuous. We further establish a martingale representation theorem for a class 
of continuous Hunt processes under certain technical conditions. In particular we 
establish the martingale representation theorem for the martingale parts of (re- 
flecting) symmetric diffusions in a bounded domain with a continuous boundary. 
Together with an approach put forward in [21], our martingale representation 
theorem is then applied to the study of initial and boundary problems for quasi- 
linear parabolic equations by using solutions to backward stochastic differential 
equations over the flltered probability space determined by reflecting diffusions 
in a bounded domain with only continuous boundary. 
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1 Introduction 



The Brownian motion, an important example of martingales, diffusions and Gaussian pro- 
cesses, possesses some remarkable properties which have been the inspiration for research in 
areas from probability, statistics to mathematical finance. 

Let B = {Bt)t>o be the Brownian motion in started with an initial distribution /i. 
The natural filtration {^l^)t>o (called the Brownian filtration, for a definition, see [5]) is 
continuous, and all martingales on the filtered probability space {Q, , P^) are contin- 

uous. More importantly the Brownian motion has the martingale representation property: 
any martingale over ^ , P^) can be expressed as an Ito integral against Brownian 
motion. 

The (predictable) martingale representation property of a family of martingales has been 
studied in a more extended setting, and several general results have been obtained. For 
example, Jacod and Yor [18] have discovered the equivalence between the martingale rep- 
resentation property and the extremal property of martingale measures. Jacod [19] also 
obtained the martingale representation property in terms of the uniqueness of some mar- 
tingale problems, and further present criteria in terms of predictable characteristics. These 
results have greatly illuminated the subject matter. When applied to specific situations, 
further work and indeed hard estimates are often required. As a matter of fact, we still have 
very limited examples of martingales and filtrations which possess martingale representation 
property (see, e.g. H], [19], [27], [33]). 

The renewed interest in recent years in the martingale presentation property has been 
motivated not only by its own right, but also by its important applications in the mathe- 
matical finance ([13], [SI]), backward stochastic differential equations ([T], [2H]) and their 
applications in some non- linear partial differential equations, see also for example [5], [B], 
[20], [21], [29], [32] and the reference therein. 

An intimate question is the continuity of natural filtrations generated by semimartingales 
and diffusion processes. A great knowledge about them has been obtained in the past. For 
example, a complete characterization of the natural filtrations generated by simple jump 
processes and Levy processes in terms of their sample paths is known. Much information 
has been obtained for a class of Markov processes. We know, from the fundamental work by 
Blumenthal, Chung, Dynkin, Getoor, Hunt, Meyer etc., (see for example [5], [9], [12], [26], 
in particular Hunt [14j), that the natural filtration of a Feller process is right continuous 
and quasi-left continuous. On the other hand, to the best knowledge of the present authors, 
there are no general conditions in literature to guarantee the martingales over the natural 
filtration of a Markov process to be continuous. A reasonable conjecture is that the natural 
filtration of a diffusion process (a continuous strong Markov process) should be continuous, 
so are the martingales over the natural filtration. Such a result is plausible but remains to 
prove. 

In this paper, we show that all martingales over the natural filtration of a continuous 
Hunt process are continuous. Our proof follows a key idea originated from Blumenthal [1], 
formulated carefully in Meyer [26] . 
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The main result of the present article is a martingale representation theorem for a class 
of continuous Hunt processes which satisfy a technical condition called the Fukushima rep- 
resentation property. 

As a consequence of our main result, we establish the martingale representation theorem 
for symmetric diffusion processes on a domain, with Dirichlet or Neumann boundary con- 
dition. More precisely, let D C M'^ be a bounded domain with a continuous boundary dD. 
Consider the symmetric diffusion in D with a formal infinitesimal generator 



subject to the Dirichlet or Neumann boundary condition (for precise meaning, see §4 and 
§5 below), where (a*-') are only Borel measurable and satisfies the uniform elliptic condition. 
To make it clear, let 



and the Dirichlet space ^ = Hq{D) or H^{D) depending on the Dirichlet or Neumann 
boundary condition. Here we use the same letter ^ to denote the filtration as well as the 
Dirichlet space: we hope it should be clear from the context which one =^ stands for. 

are the natural filtrations generated by the symmetric diffusion {Xt)t>Q- It happens in 
this case that the coordinate functions u^{x) = belong to the local Dirichlet space ^loc, 
and Xf = {Xl , • • • , X^) has the Fukushima's decomposition 



for all a; G -D (or D in the Neumann boundary condition case) except for a zero capacity set 
with respect to the Dirichlet form [S',^), where M^, ■ ■ ■ , are continuous martingales 
additive functionals and A^, • • • , A'^ are continuous additive functionals with zero energy. 
The following martingale representation theorem follows from our main result. 

Theorem 1.1 Under the above assumptions, for any initial distribution ji which has no 
charge on capacity zero sets, the family {M^ , ■ ■ ■ , M"') of martingales over (fi, P^) 
has the martingale representation property: for any square integrable martingale N — {Nt)t>o 
over (fi, ^l^, P^) there are unique {^^) -predictable processes , • ■ ■ , F'^ such that 






Xl-Xl^Ml+Al P--a.e. j = 1, • • • , 



d 
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Theorem 11.11 may be applied to the symmetric diffusions in with Dirichlet form ^) 
where ^ = H^{M.'^). This special case has been proved in [33] and [1]. 

In [33], Zheng has pointed out that the martingale part of symmetric diffusion {Xt)t>o in 

with infinitesimal generator being a uniform second order elliptic operator in divergence 
form has the martingale predictable representation property and described a proof based 
on the results on the Dirichlet process p{t,Xt) obtained in Lyons and Zheng [22| and |23j . 
where p{t, x) is the probability density function of Xt under the stationary distribution. 
More precisely, Lyons and Zheng [23j have extended Fukushima's representation theorem for 
martingale additive functionals to a class of processes which has a form f{t, Xt), where / has 
finite space-time energy. Their results in particular yield that p{t, Xt) is a Dirichlet process 
in the sense of Follmer [T5|, and its martingale part can be expressed as an Ito integral 
against {M^, ■ ■ ■ , M'^), which, together with the Markov property, allows to show that for 
^ = fi{Xt-^)---fn{Xt„) the conditional expectation E{^\^t) is again a Dirichlet process 
which can be expressed as an Ito's integral against (M^, ■ ■ ■ , M'^), where the expectation is 
taken against the stationary distribution P™(-) = J^^ P^[-)dx. A routine procedure based 
on the Doob's maximal inequality allows to prove the martingale representation theorem 
for the symmetric diffusion in with the generator L an elliptic operator of second order. 
Apparently not knowing the work [33|, in an independent work Bally, Pardoux and Stoica 
|TJ, among other things, a detailed proof has been provided. 

The technical difficulty with the proof described above lies in the fact that even for a 
smooth function /, f{Xt) may not be a semimartingale, so that Ito's calculus can not be ap- 
plied. Instead of considering random variables with product form such as = fi{Xt^) ■ ■ ■ fn{Xt„) 
which linearly span a vector space dense in Lp{Q, P^), we utilize a linear vector ^ spanned 
by those random variables which have a product form ^ = Ci' ' ' with 

POO 

= / e~'^^'f,{Xt)dt 
Jo 

where fj are bounded Borel measurable functions, and aj > for j = 1, ■ ■ ■ ,n. According 
to Meyer ^ is dense in L^{Q, ^ , P^^). The important feature is that U'^f{Xt) is a 
semimartingale for any a > and a bounded Borel function /, where is the resolvent 
of the transition semigroup. Moreover, in the symmetric case, always belongs to the 
Dirichlet space (when / is bounded and square integrable) and thus Fukushima's representa- 
tion theorem for martingale additive functionals can be applied to extend the representation 
to any martingales. 

The martingale representation theorem for the symmetric diffusion process {Xt)t>o in a 
domain D allows us to study the following type of backward stochastic differential equation 
(BSDE) 

d 

dYt = -fit, Yt, Zt)dt + J2 ZidMi 

i=i 

with a terminal condition that 1^ = 1^^ L'^i^-, P^)) and thus gives a probability rep- 
resentation for weak solutions to the initial and boundary value problem for non-linear 
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parabolic equations. The existence and uniqueness of solutions to the BSDE follows from 
exactly the same approach as for the Brownian motion case, which is the pioneering work in 
BSDE done by Pardoux and Peng in [28]. We however describe an approach put forward in 
pi] . which allows us to devise an alternative probability representation for the initial and 
boundary problem of the corresponding semi-linear parabolic equation 



subject to ■^iQjjU^t, ■) = in a bounded domain with only continuous boundary. 

The paper is organized as follows. In Section 2, we develop further an idea put forward 
in Meyer [26], and show both the natural filtrations and martingales over the natural filtra- 
tions for a continuous Hunt process are continuous. This result is hardly new but it seems 
not appear in the literature yet. In order to prove this result, we devise an important while 
elementary formula for which may be considered as a refined version of a classical 

formula devised firstly by Hunt and Blumenthal for potentials and multiple potential case by 
Meyer. In Section 3, we establish the main result of the paper: a martingale representation 
theorem for a continuous Hunt process under technical assumptions called the Fukushima 
representation property, and give some examples in which our result may apply. In Section 
4, we outline the existence and uniqueness of solutions to backward stochastic differential 
equations over the natural filtered probability space over a reflecting symmetric diffusion in 
a bounded domain with non-smooth diffusion coefficients and non-smooth boundary, and 
finally we apply the theory of BSDE to the study of the initial and (Neumann) bound- 
ary problem of a non-linear parabolic equation in a bounded domain with only continuous 
boundary. We believe these results are new even for reflecting Brownian motion in a domain 
with non-smooth boundary. 

2 Martingales over the filtrations of continuous Hunt 
processes 

Consider a Markov process {Q, ^j^, Xt, 9t, P^) in a state space E' = EU {d}, where E is 
a locally compact separable metric space E, with transition probability function {Fj(a;, ■) : 
t > 0}, i.e.. 



JE 

In (12.11) . E^ on the left-hand side stands for the (conditional) expectation with respect to 
the probability measure P^, and the right-hand side may be abbreviated as Ptf{Xs) where 



JE 

which is well defined for a bounded or non-negative Borel measurable function /. The family 
of kernels {Pt)t>o is called the transition semigroup associated with the Markov process 
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{Xt)t>Q- Without specification, (^°)t>o is the filtration generated by {Xt)t>o, that is = 
a{Xs : s<t} for t > and = a{Xs : s > 0}. 

For a cr-finite measure // on {E,^{E)) (where ^{M) always represents the Borel a- 
algebra on a topological space M) 



defines a measure on ^^). If /i is a probability, then {Xt)t>Q is Markovian under P'^ with 
transition semigroup {Pt)t>Q and initial distribution jj,, in the sense that 



and P'^IXq e A} = ii{A) for any A e ^{E), where E^^ is the (conditional) expectation 
against the probability P'^. 

Denote by ^{E) the space of all probability measures on {E,j3§{E)). If /i G ^{E), 
denotes the completion of under P^, and is the smallest a-algebra containing 
^^ and all sets in with zero probabihty. (^/^)t>o is called the natural filtration of the 
Markov process (Xt)t>o with initial distribution yu. Let 



If (^t)t>o is a filtration, i.e. an increasing family of cr-algebras on a common sample space, 
then = '^s>t^s for i > and %- = a\^s : s < t} for t > 0. The filtration is called 
right (resp. left) continuous if = for all t > (resp. '^t- = for allt > 0). The 
sample function properties of a Markov process (Xt)t>o and the continuity properties of its 
natural filtration had been studied by Blumcnthal, Dynkin, Getoor, Hunt, Meyer etc. The 
fundamental results have been established via the regularity of the transition probability 
function {Pt(x, •) : t > 0}. Their work achieved the climax for Markov processes with Feller 
transition semigroups. 

As matter of fact, the continuity of the filtration (i^/') (or (#*)) does not follow that of 
sample function {X^i~^q. For example, a right continuous Markov process docs not necessarily 
lead to the right continuity of its natural filtration i^^i ) (or {^^t))- The same claim applies 
to the left continuity. In fact, the regularity of natural filtrations is much to do with the 
nature of the Markov property, such as strong Markov property. 

Let Coo(P) (resp. Cq(E)) denote the space of all continuous functions f on E which 
vanish at infinity d, i.e. lim^^a/(^) = (resp. with compact support). Recall that a 
transition semigroup {Pt)t>o on {E,^{E)) is Feller, if for each t > 0, Pt preserves Coo{E) 
and limtj^o Ptf{x) = f{x) for each x E E and / e Coo{E). 





and (^ 
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For a given Feller semigroup {Pt)t>o on {E,^{E)), there is a Markov process 

with the Feller transition semigroup {Pt)t>o such that the sample function f — > is right con- 
tinuous on [0, oo) with left hand limits on (0, oo). In this case, we call {Q, Xt, 6t, P^) 
a Feller process on E. 

For a Feller process (fi, X^, 6*4, P"^), the natural filtration {<^l^)t>Q for any E 

l3^{E) and as well as {^t)t>o are right continuous. {Xt)t>o and {^t)t>o are also quasi-left 
continuous, that is, if is an increasing family of (t^/^)-stopping times, and T„ | T, then 
\\mn^ooXT„ = Xt on {T < 00} and = a{^^^ : n G N}. Therefore accessible {^t)- 
stopping times are predictable. An (^/')-stopping time T is totally inaccessible if and only 
if P'^{T < 00} > and Xt 7^ Xt- on {T < 00} P^-a.e. Similarly, T is accessible if and 
only if Xt = Xt- P^-a.e. on {T < 00}. Hence X has only inaccessible jump times. 

What we are mainly concerned in this article is Hunt processes. Hunt processes are right 
continuous, strong Markov processes which are quasi-left continuous. These processes are 
defined in terms of sample functions, rather than transition semigroups, see |5] and ^ for 
details. It is well-known that Feller processes are stereotype of Hunt processes or the later 
is an abstraction of the former. 

We are interested in the martingales over the filtered probability space {Q, t^^, P^), 
and we are going to show that, if {Q, Xt, Ot, P^) is a Hunt process which has contin- 

uous sample function, then any martingale on this filtered probability space is continuous, 
a result one could expect for the natural filtration of a diffusion process. Indeed, this result 
was proved more or less by Meyer in his Lecture Notes in Mathematics 26, "Processus de 
Markov" . Meyer himself credited his proof to Blumenthal and Getoor, more precisely a cal- 
culation done by Blumenthal However it is surprising that the full computation, which 
yields more information about martingales over the natural filtration of a Hunt process, was 
not reproduced either in the new edition of Meyer's "Probabilites et Potentiels" or Chung's 
"Lectures from Markov Processes to Brownian Motion", although it was mentioned in [TU] 
where Chung and Walsh gave an alternative proof of Meyer's predictability result, so that 
Blumenthal's computation is no longer needed. However, it is fortunate that Blumenthal's 
calculation indeed leads to a proof of a martingale representation theorem we are going to 
establish for certain Hunt processes, see section §3. 

Let us first describe an elementary calculation, originally according to Meyer [SB] due to 
Blumenthal. Let 

X = {n,^\^;',Xt,et,pn 

be a Hunt process in a state space E' = EU{d} with the transition semigroup {Pt)t>o, where 
d plays a role of cemetery. Let : a > 0} be the resolvent of the transition semigroup 
{Pt)t>o- 

POO 

U''{x,A) = / e''''Ptix,A)dt 
Jo 
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and (f/")a>o the corresponding resolvent (operators), i.e. 

f/-/(x) = / fiz)U''ix,dz) 
Je 

POO 

= / e-'''Ptf{x)dt 
Jo 

for bounded or nonnegative Borel measurable function / on E. To save words, we use 
^h{E) to denote the algebra of all bounded Borel measurable functions on E. Obviously, 
CUE) C ^b{E). 

Let K{E) C j3§f,{E) be a vector space which generates the Borel cx-algebra J3§{E). Let 
C L^{fl, P'^) (for any initial distribution /x) be the vector space spanned by all ^ = 
^i - ■ ■ for some n G N, 

/•oo 

= / e-^^'f,{Xt)dt 
Jo 

where aj are positive numbers, fj G K{E), j = 1, - ■ ■ ,n. Meyer [26] proved that ^ is dense 
in L^{Q, P^) for a Hunt process. Since this density result will play a crucial role in what 
follows, we include Meyer's a proof for completeness and for the convenience of the reader. 
The key observation in the proof is the following result from real analysis. 

Lemma 2.1 Let T > 0. Let IK denote the vector space spanned by all functions ea{t) = e~°*, 
where a > 0, then IK is dense in C[0,T] equipped with the uniform norm. 

The lemma follows from Stone- Weierstrass' theorem. 

Lemma 2.2 (P. A. Meyer) For any initial distribution ^ and p G [l,oo), ^ is dense in 
LP(fi,^^,P^). 

Proof. First, by utilizing Doob's martingale convergence theorem, it is easy to show 
that the collection s^/ of all random variables which have the following form 

9i{Xt,) ■ --gniXtJ, 

where n G N, < ti < ■ ■ ■ < t„ < oo and G K{E), is dense in LP{n, ^f", P^). Let M' be 
the linear space spanned by all = ?7i ■ ■ ■ rjn, where 

POO 

Vj = / 9j{Xt)^j{t)dt, 
Jo 

where Qj G K{E) and ipj G C[0, 00) with compact supports. According to the previous 
lemma, for every e > we may choose ipj G K such that 

\(Pj{t) - ipj{t)\ < ee"^* for alH > 
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for some A > 0. Let 

/■oo 

= / 9jiXt)ilj,{t)dt. 
Jo 

Then ^ = ^i - ■ ■ ^ and 

\Vj-^j\ < ^ll^illoo^ 

where || ■ ||oo is the supermum norm. It follows that 

i?ie-er<-max||^,ii^%^ 

and thus ^ belongs to the closure of Finally it is clear that any element 

giiXtJ ■ ■ ■ g^iXtJ = lim Vi---Vn 

fe— >oo 

where 



ri] = / g,iX,)ip^it)dt 
Jo 

and ip^ has compact support and ip^ 6t^ weakly. We thus have completed the proof. ■ 

Let be any fixed initial distribution. If / is a bounded Borel measurable function on 
E and a > 0, then 

poo 

e= / e-'^'f{Xt)dteL\n,^^,P^). 
Jo 

Consider the martingale Mt = E'^ {^|^/^} where t > 0. According to an elementary compu- 
tation in the theory of Markov processes, 

pt f POO 

= j e-''J{Xs)ds + E>'lj^ e-'''f{Xs)ds\^ 

pt POO 

= / e-''JiXs)ds + e-"* / e-''-'PsfiXt)ds 
Jo Jo 



It is known that if X = {Xt)t>o is a Hunt process, then for any a > and bounded Borel 
measurable function /, If^f is finely continuous, i.e., t U^f{Xt) is right continuous. 
Moreover if X is a continuous Hunt process, it follows from a result proved by Meyer that 
t — i> U"f{Xt) is continuous, and therefore, the martingale Mt = E^{^\^t^ is continuous. 
We record Meyer's result as a lemma here. This result was proved in [26] for Hunt processes 
(see T15 THEOREME, page 89, ^26]). A simpler proof for Feller processes may be found on 
page 168, [H]. 
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Lemma 2.3 (Meyer) Let {Xt)t>o be a Hunt process, f e ^b{E), a > and h — U°- f he a 
potential. Then 

h{Xt-) = h{X)t- Vt > P^'-a.e. 

for any initial distribution ^. 

P. A. Meyer pointed out that the previous computation can be carried out equally for 
random variables on (fi, which have a product form ^ = ^i - ■ ■ in where each — 
e^°'i^fj{Xs)ds. Let 7r„ denote the permutation group of {!,••• ,n}. 

Lemma 2.4 (Blumenthal and Meyer) Let ^ = Ci---^n where = e~'^i^ fj{Xs)ds, 

aj > and fj e ^b{E), and Mt = E^' {i\^t}. Then 

Mt = Yl E (tlf e-"-7..(^.)^^ ) • F^n,-,J.-,Jn)m (2.2) 

k=0 (jl,- ,jk,- ,jn)e7rn \i=i ''^ J 

where 

F(Ju-,j.-,a^) = I n e-'^nsf.^x^yi,^ | . (2.3) 

Proof. The task is to calculate the conditional expectation Mt = {i\^t}- The idea 
is very simple: spliting each into 

pt l-CO 

ij^ e-"^V,(^.)c^s + e-"^* / e-'^i'f^{Xsoet)ds 
Jo Jo 

so that 



fc=0 (il,- Jfe,- Jn)67r{l,-,n} \i=l-^° 

( Ji- poo \ 

n ^"^'''l e-'^n^f,,{X,oet)ds\. 



By using the Markov property one thus obtains 
Mt = E^{i\^t} 



= E E n/ 

•E^ I || n e-'^n^ J\-''n^fj,{X, o 0t)ds^ 

= E E (ufe--^hiXs)ds) 



■E 



•Xt 



^ poo \ I 

n I e--^''U{X,)ds\\. (2.4) 
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Our only contribution in this aspect is the following formula, which allows to prove 
not only that all martingales over the natural filtration of a continuous Hunt process are 
continuous, but also a martingale representation theorem in the next section. 

Lemma 2.5 Let aj be positive numbers and fj e ^b{E) for j — 1, ■ ■ ■ , k. Consider 
F{x) = [■■■[ e-^U-^^^ [ h{z^)---h{zk)PsAx,dz^) 

J Jo<si<-<Sk«xi J E®'' 

■Ps^-siizi, dz2) ■ ■ ■ Psi^-sk-ii^k-i, dzk)dsi ■ ■ ■ dsk. 

Then 

F = (/i(t/°2+-+°'=/2 ■ ■ ■ (U^^fk) ■■■) . (2.5) 

Proof. To sec why it is true, we begin with the case that A; = 1. In this case F — 
e-'^'Psfds = mf. If A; = 2, then 



F ^ e-'''''-'e-^^'^Ps,{hPs,-sj2)ds,ds2 

J J 0<Si<S2<OO 
noo poo 

= / / e-''^'e-''^'Pt{hPs.tf2)dsdt 
Jo Jt 

= e-'^^'Pt {^j^ e-'^^' hPs-tf2d^ dt 



OO / pOG 

— Oi\t —OL2t 



^^a,t^-a,tp^ (^/l j e-''''Psf2ds ] dt 



OO 
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and by an induction argument, for a general case. Indeed, ii k > 2, then 



F{x) = I ■ ■ ■ I e-^^"=i"^'^e 

J ^0<si<---<Sfc_|_i<oo 



xP,j [x, dzi)Ps^-si{zi, dz2) ■ ■ ■ Psi^-sk-ii, Zk-i, dzk)dsi ■ ■ ■ dskdsk+i 



0<si<---<S(t<oo 
oo 

xPs^{x, dzi)Ps^-si{zi,dz2) ■ --Ps^-sk-ii, Zk-i,dzk)dsi ■■■dsk 



0<si<---<Sfe<oo Je^fc 



X (^/fe(^fe)e-"'=+i^'=^ e-''''+^'Ptfk+i{zk)dt^ 

xPs^{x, dzi)Ps2-siizi, dz2) ■ ■ ■ Psk-sk-ii,Zk-i,dzk)dsi ■■■dsk 



e 

0<si<---<Sfc<cx) 



X / fM)---h{zk)U->'+^h+iizk)dt 

and the formula follows the induction assumption. ■ 
Lemma 2.6 Let /i, ■ ■ ■ , /fc G ^h{E), aj positive numbers, and 

{I /"OO 

Then 

F = J2 U'^''^-^""' (/.il^"^'^^-^"^-^/.. • • • {U'-'^fjk) ■■■) (2-6) 
O'l,— ,jk}^^k 

where tt^ zs i/ie permutation group of {!,■ ■ ■ ,k}. 
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Proof. We have 



f POO POO ^ 

POO POO 

Jo Jo 



/■■■/ e-"i«i...e-°fe«. 

{ju- J0<s,^<-<s,^<oo 

xE-{MXj---h{Xj}ds,---dsk 

/ ■ ■ ■ / e""^!''^ ■ ■ ■ e'^'k^k 

r- 1^ >/ 0<si< - <si-<oo 

{h,-Men,y J0<s,<-<s,<coJe^>' 

xPs,{x, dzi)Ps^-siizi, dz2) ■ ■ ■ Ps^^-sk_,izk-i,dzk)dsi ■■■dsk 



and (12.61) follows from Lemma 
From now on, we assume that 



X = {n,^',^}>,Xt,9t,p-) 

is a continuous Hunt process in E' = E U {d} with the transition semigroup {Pt)t>o- In 
other words, it is a Hunt process with continuous sample paths. Therefore, {Xt)t>o is a 
diffusion process in E, i.e. {Xt)t>o possesses the strong Markov property with continuous 
sample function. Under our assumptions, any finite (.^/')-stopping time is accessible and 
thus predictable, and therefore = In particular, {^t) is left continuous, and thus 

the filtration (^Z') is continuous for any initial distribution /i. 

Since any martingale on t^^, , P'^) has a right continuous modification, by a mar- 
tingale we always mean a martingale with right continuous sample function. 

Lemma 2.7 Suppose ^ = ■ ■ ■ where each has the following form 

POO 

= / e-''^J,iX,)ds 



where aj > and fj G J!^b{E)- Let Mt = E'^ {i\'^t}- Then {Mt)t>o is a bounded continuous 
martingale on {Q, , P^) . 

Proof. According to Lemma [231 we need only to show that for function of the following 
type 



F(x) = E'^ 



" POO \ j 

n e""^'* / e-'^n^f,,iX,)ds] 
:=k+i / J 
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t F{Xt) is continuous. By Lemma [231 F is an a-potential, so that it is finely continuous, 
and together with Lemma [2. 3[ it imphes that t — > F{Xt) is continuous, which completes the 
proof. ■ 

We now state the main result of this section. For simplicity, a square integrable martingale 
{Mt)t>o over (fi,^,^t,P) means Mt = E{^\^t) with ^ G L'^iVt,^ ,P). This is equivalent 
to say supj^o ^[^t] < 

Theorem 2.8 Let 

he a continuous Hunt process in E, and fi E 0^{E). If ^ e L2(1], P^), then the martin- 
gale Mt = E'^ {^1^/*} is continuous, that is, square-integrable martingales on {Q, ^1^, P^) 
are continuous. Therefore local martingales on {Q, , P^) are continuous. 

Proof. We can choose a sequence ^„ G ^ such that ^„ — > ,^ in L^. Doob's max- 
imal inequality implies that, if necessary by considering a subsequence, the martingales 
{E^{^n\'^t) : ^ > 0} converges (almost surely at least along a subsequence) to {E^^^l^l") : 
t > 0} uniformly on any finite interval of t > 0. It is shown in Lemma 12.71 that for 
each n, the martingale E'^{^n\'^t) is continuous and thus the square integrable martingale 
{E^'{^\^l') : t > 0} must be continuous. 

By the localization technique, it follows thus that local martingales on ^l", P^) 

are continuous. ■ 



3 Martingale representation for continuous Hunt pro- 
cess 

In this section we assume that 

x = {n,^^,^^,Xt,et,pn 

is a continuous Hunt process in the state space E' = E U {d} with transition semigroup 
{Pt{x,dy) : t > 0}, where i? is a locally compact separable metric space. Let /i G ^{E) be 
an initial distribution. 

If a > and / G ^b{E) then M"''^ denotes the continuous martingale 

Mt^ = e-''J{Xs)ds\^r^ . 

Recall that, if u is an a-potential, i.e., u = U^f where / G j3§b{E), then u{Xt) — u{Xq) is a 
continuous semimartingale on {^'^ , , P'^). and possesses Doob-Meyer's decomposition 

u{Xt) - u{Xo) = + 
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where ^ ^ 

mJ"1 = / e'^'dM^'^, 4"' = / Lu{X,)ds 
Jo Jo 

and Lu = au — f. 

We make the following assumptions on the continuous Hunt process X started with an 
initial distribution G ^{E), and we call these assumptions the Fukushima representation 
property. 

Assumptions. There is an algebra (a vector space which is closed under the multiplication 
of functions) K{E) C ^b{E) which generates the Borel cr-algebra ^{E) and is invariant 
under [/" for a > 0, and there are finite many continuous martingales M^,- ■ ■ , M'^ over 
{Q, , P^) such that the following conditions are satisfied: 

(1) For any potential u = U"f where a > and / G K{E), the martingale part Af'^l 
of the semimartingale u{Xt) — u{Xo) has the martingale representation in terms of 
(M^, ■ ■ • , M'^), that is, there are predictable processes Fi, ■ ■ ■ , Faon {Q, such 
that 

d i-t 

M'f = V / FidMi P^-a.e. (3.1) 

(2) {{M^,M^)t) is strictly positive definite. 

The first assumption means that the martingale M'"! with u being a potential may be 
represented. The second condition ensures that the representation (13. ip is unique. 

The Fukushima representation property is mainly an abstraction of the chain role for the 
martingale part of u{Xt). Indeed, if Xt = {X}, ■ ■ ■ , X^) is a d- dimensional Brownian motion 
and u is an a-potential with a > 0, then u is smooth and by Ito's formula 

^ pi o pt "1 

u{Xt)-u{Xo) = J2j^ ^{Xs)dXi + -Au{Xs)ds 

so that 

One can easily see that the Brownian motion satisfies the Fukushima representation property. 

Theorem 3.1 (Martingale representation) Let^E^{E). Suppose that the Fukushima 
representation property holds for X with a finite set of martingales {M^, ■ ■ ■ , M'^). For any 
square-integrable martingale N = {Nt)t>o on {Q,^^,^I^,P^), there are unique predictable 
processes (F/) such that 

Nt-No = ^ I FidMi P^-a.e. 
i=l ^0 
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Proof. The uniqueness follows from condition 2) in the Fukushima representation. We 
prove the existence. Take ^ e L'^{Q, P^) such that = E''{^\^l'}. Since ^ is dense in 
L'^{Q,^'^, P^), so we first prove the martingale representation for ^ G By the linearity, 
we only need to consider the case that C = ^i - ■ ■ Cn where = e'""^^ fj{Xs)ds for aj > 
and fj G K{E). In this case, according to \2A\ Lemma [2.51 and Lemma [2.61 

m 

where the sum is a finite one, and for each m, Z"^ = Zt has the following form 

(the superscript m will be dropped if no confusion may arise), where 

Vt = ll [ e-^^'g,{X,)ds 
i=i 

and 

- [■■■[ e-^Uf^^^^ [ h{zi)---hk{zk)PsAx,dz,) 

J Jo<si<--<Sk<oo J E®^ 

xPs2-si{zi,dz2) ■ ■ ■ Ps^-sj,_^{zk-i,dzk)dsi ■■■dsk 

for some k' and k, (3i > and functions gi, hj are bounded and continuous. According to 
Lemma 12.51 

u = f/ft+-+A (/ii(?7^2+-+^^/i2 ■ ■ • (U^^hk) ■■■). 
In particular, u is again a potential which has a form u = If^g for 

g = hi{U^^+-+^^h2 ■ ■ ■ {W^'hk) ■■■)e K{E) 

and a = Pi + ■ ■ ■ + l3k- Hence u{Xt) is a continuous semimartingale with decomposition 

where A^"! is continuous with finite variation, and due to the Fukushima representation 
property 



uix] 



7 = 1 



for some predictable processes GK In particular, each is a continuous semimartingale. 
Since, by Theorem 12.81 is a continuous martingale, so that 

Nt = the continuous martingale part of V^v!^{Xt). 

m 
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Therefore we are interested in the martingale part of Zt — Vtu{X-^. Since F is a finite 
variation process, so according to Ito's formula 



Jo 

^ Zo+ [ u{X,)dVs+ I VsdAY^+ f Kc/m'"' 
Jo Jo Jo 

Jo Jo j^^i Jo 

so that the martingale part of Zt is 

J2 fvs-G\dMl. 



Therefore 



which shows the martingale representation. 

Suppose now ^ G LF'iVL, , P'^). Choose a sequence ^„ G such that ^„ — >• ^ in 
L'^{n,^^',P^'). Let iV^^"' = and Nt = E^'{i\^^). According to Doob's maximal 

inequality, if necessary by passing to a subsequence, we can assume that N^^^ converges to 
Nt uniformly on any finite interval. Nt^^ has the martingale representation 

d rt 



.7=1 



so that 

An) ^j[m) ^j{n) ^j[m) 



J2 f\F{n)l- F{m)l){F{n)i- F{m)i)d{M\M^),. 



Since {{M\M^)t) is positive, it follows that (F(n)^,--- ,F{nY) converges to predictable 
processes (F^, • • • , F'^) under the norm 

CO ^ d Y fN 

ll(^\---,^')ll = E^E^' / KFid{M\M^)s . (3.2) 
Ar=i i,j=i L^o -I 

Then 

d 



Nt-No = j2 f ndMi 

7 = 1 ^0 



3= 
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This theorem claims that as long as every martingale of resolvent type is represent able, 
so is any martingale. When is the Fukushima representation property satisfied? There are 
many examples. In the remain of this section, we shall give three interesting examples in 
symmetric situation. 

Brownian motion with any initial distribution is certainly an example. Indeed, for Brow- 
nian motion in Mf^, we may choose K{E) = C^{R'^) (the space of smooth functions which 
vanish at infinity), then for / G K{E), If^f is smooth, and (13.11) follows from Ito's formula 
applying to U^f. Theorem 13.11 gives a new proof for classical martingale representation 
theorem. 

The second example is the reflecting Brownian motion. As Example 1.6.1 in [T7], we 
consider Dirichlet form (|D,i/^(D)) on L'^{D) where D is the classical Dirichlet integral 
and D is a bounded domain on R''. We further assume that any x G dD has a neighborhood 
U such that 

DnU = {{Xi) eR'':Xd> F(xi, ■ • ■ , Xrf_i)} n U 

for some continuous function F. Then C^{D) (the space of restriction to D of functions in 
C^(R'^)) is dense in H^{D) (see [22] for details), i.e., (|D, H^{D)) is a regular Dirichlet form 
on LF'{D). The corresponding continuous Hunt process X = {Xt, P^) is called the reflecting 
Brownian motion. For x = (x*) G M.'^, we use Ui{x) = x'^ , 1 < i < d, to denote the coordinate 
functions. Then Ui E ^ and we denote by M* = M''*'! the martingale part in Fukushima's 
decomposition. It can be seen from Corollary 5.6.2 [T7] that for any u G C^{D), 

Mf^ = V / ^(Xs)dMi, P"-a.s. for q.e. x eD, 
j^Jo dxi 

where q.e. means 'quasi-every where', i.e., except a set of zero-capacity. Then a routine ap- 
proximation procedure shows that for any u G H^{D), there exist Borel measurable functions 
{fi : 1 < i < d} on D such that 

^ = / fi{Xs)dMl P^-a.s. for q.e. x G R'^. 
i=l ^0 

Therefore the reflecting Brownian motion has Fukushima representation property, by choos- 
ing K{D) to be the space of bounded measurable functions and any initial distribution ^ 
charging no set of zero capacity, i.e., a smooth distribution, because an exceptional set exists 
in above representation as is always when the process is constructed through a Dirichlet 
form. If the boundary is Lipschitz, then the transition function has density ([2j) and in 
this case, the exceptional set may be erased. Notice that under the current condition, the 
reflecting Brownian motion X itself is not necessarily a semimartingale. The readers who 
are interested may refer to [2], [Tj and [8] about when a reflecting BM is a semimartingale 
and the corresponding Skorohod decomposition. It should be pointed out that, although the 
martingale part of the reflected Brownian motion is a Brownian motion, but the martingale 
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representation property does not follow from the classical representation property for Brow- 
nian motion. The reason is that, as long as the boundary is not sufficiently smooth, the 
natural filtration {^t)t>o is much bigger in general than the natural filtration generated by 
the martingale part (M^, ■ ■ ■ , M'^) of X. 

Another example our main result may apply is symmetric diffusions in a domain killed 
at boundary. Actually Theorem 6.2.2 in [T7j tells us that every continuous symmetric Hunt 
process with a smooth core enjoys the Fukushima representation property. More precisely 
let D be a domain of with continuous boundary dD and m a Radon measure on D. 
Let X be a continuous Hunt process which is symmetric with respect to m and ^) the 
associated Dirichlet form on L'^[D,m), which has Cl{D) as a core. For x = (x*) G M'^, we 
use Ui{x) = x"^ , 1 < i < d, to denote the coordinate functions. Then Ui G ^loc and we denote 
by M* = M'"'] the martingale part in Fukushima's decomposition. Let 

the smooth measure associated with CAF {M\ M^). Then <f is expressed as 

As asserted in Theorem 6.2.2 [17], for any initial smooth distribution fi (i.e. a probability 
on (D, ^{D)) having no charge on capacity zero sets) and u G the martingale part Afl^l 
in Fukushima's decomposition of u may be represented as 

^1"'=E f fi{Xs)dMi P'^-a.e. 

where /i, ■ ■ ■ , /d G SS{D). If we take K{E) = L'^{E, m) n ^b{D), X satisfies the Fukushima 
representation property. In these examples, {M*} are the martingales corresponding to 
coordinate functions so we call them coordinate martingales. 

To have the uniqueness, some kind of non-degenerateness is needed. We say that X 
is non-degenerate if the condition (2) in Fukushima representation property is satisfied: 
{{M\M^))i<ij<d is positive. 

Corollary 3.2 Assume that X is either the reflecting Brownian motion on a bounded do- 
main or a non-degenerate symmetric Hunt diffusion on a domain D (ZW^ as stated above. 
Then the Fukushima representation property is satisfied and therefore the martingale rep- 
resentation holds in the sense of Theorem \3.1\ with coordinate martingales and for a given 
initial distribution fi charging no sets of zero capacity. 

From this result, we may recover the martingale representation established in [Tj and [33], 
where X is a diffusion process corresponding to non-degenerate symmetric elliptic operator 
on W^. 

Without essential difference, the conclusion holds also for reflecting diffusions on such 
domain with generator being a symmetric uniformly elliptic differential operator of second 
order as introduced in the beginning of next section. 
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4 Backward stochastic differential equations 



In this section we consider backward stochastic differential equations which can be used 
to provide probabihty representations for weak solutions of the initial and boundary value 
problem of a quasi-linear parabolic equation. 

Let D C M"^ be a bounded domain with a continuous boundary dD^ D = D U dD the 
closure of D. Let 

L = i V —a'^ — 
2 ^ dx^ dx' 

be an elliptic differential operator of second order, where a = (a^^) is a positive-definite, 
symmetric, matrix-valued function on D, a = (a*-') is Borel measurable, and satisfies the 
elliptic condition: 

d 

for all x G -D for some constant A > 0. Consider the Dirichlet form [S",^) on L'^{D,dx), 
where 

and ^ = H\D). 

Let Q he a. space of all continuous paths in D, {Xt)t>o the coordinate process on f2, 
= (7{Xs : s > 0}, = (j{Xs '■ s < t} for each t > 0, and {Ot)t>o shift operators on Q. 

Let 

be the canonical realization of the symmetric diffusion process in the state space D associated 
with the Dirichlet space {S',^), which is called a reflecting symmetric diffusion in D. 

The coordinate functions Uj{x) = x^ {j = 1, ■ ■ ■ ,d) belong to the local Dirichlet space 
^loc, so that 

Xi - X', = Ml + A^, P^-a.e. j = 1, ■ ■ ■ , (4.2) 

for all a; G Z) except for a capacity zero set, where = M^^^^ etc. 

Let ^i{D) denote the space of all probability /i G ^{D) which has no charge on zero 
capacity sets (with respect to the Dirichlet form {<S',H^{D)) defined by (14.11) . According 
to Theorem 13.11 for any initial distribution G S^i{D), the family of martingales {M^ : 
j = I,-- - ,d} over [Q, , P^^) has the martingale representation property: for any 

square-integrable martingale = {Nt)t>Q on {Q, , , P^^), there are unique predictable 
processes (F/) such that 

t 

iVt - iVn = 2^ / FldMl F^-a.e. 



1=1 
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Let us work with a fixed smooth initial distribution /i G S^i{D) and the filtered probability 
space {n,^^',^|',P^'). 

Consider the following backward stochastic differential equation 

d 

dY; = -fit, Y„ Z,)dt + Zl'dMl = C (4.3) 

i = 1, - ■ ■ ,d', where T > 0, G L^(fi, , P^) are given terminal values, and /* are Lipschitz 
functions: there is a constant Ci > 

\r{t,y,z)\<C,{l + t+\y\ + \z\) 

and 

\P(t,y,z)-^{t,y,~z)\<C^{\y-y\ + \z~~z\) 

for all t > 0, y G M*^ , -2, 5 G M.'^ One seeks for a solution pair (y, Z) which solves the 
following integral equation 

d pj' 



/ r{s,Y„Z,)ds-J2 [ Z'JdMl (4.4) 



for t G [0, T]. The integral equation (14.41) has a unique solution pair (F, Z) such that is a 
continuous semimartingale, and Z^^ are predictable processes satisfying 



pT d 

Yl o!'\Xs)Z'^Zfds < oo. 



This can be demonstrated by employing the Picard iteration for {Y. Z) as in the case of 
Brownian motion (see [2S])- Another approach, proposed in a paper by Lyons, Liang and 
Qian [2T] which applies to a general filtered probability space, may be described as follows. 
The idea is to rewrite the integral equation (14.41) into a functional differential equation for 
the variation process part V of Y. Let Y = N — V where y is a finite variation process, and 



7 = 1 



On the other hand 

Nt = E>'{i + VT\^t]- 

Since y is a continuous semimartingale, its decomposition is unique up to an initial value. 
The integral equation (14.41) leads to that 



yt = -j^ fis,Y,,Z,)ds + NT-^ 
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conditioned on and we obtain 



= -E^ y^ f{s, Ys, Z,)ds\^t^ + {Vrl^n 
= -E'^y^ r{s,Ys,Z,)ds-VT\.^ty f{s,Ys,Z,)ds. 
Therefore the integral equation (14.41) is equivalent to 

Vt-Vo= [ f{s,Y,,Z,)ds (4.5) 



where 

Yt = Y{V)t = N{V)t - Vt, NiV)t = E'^{C + VtI^} 
and Zt = Z{V)t is determined by the martingale representation theorem 

N{v)i - N{v% = J2 f ziyr^dMi. 

Equation (14.51) thus may be written as a functional equation 

Vt-Vo= [ fis,Y{V)s,Z{V)s)ds (4.6) 



where Y{V) and Z{V) are considered as functionals of V. The Picard iteration applies to 
(14. 6 p we have 

Theorem 4.1 If ^ E L^(fi, ^^jt, P^) and /* are Lipschitz continuous, then there is a unique 
pair {Y, Z) such that Y is a continuous semimartingale which solves BSDE (^7^. 



For a complete proof of Theorem 14. ![ the reader may refer to [21 



5 Non-linear parabolic equations 

We are under the same setting as in the previous section, and use the notations established 
therein. 

To motivate our approach, let us begin with the case that a is smooth, and D is bounded 
domain with a smooth boundary. 

In this case = {Xj)t>o in (14. 2 p are continuous semimartingales, thus are finite 
variation processes. For any h G C^'^([0, oo)xD) satisfying the Neumann boundary condition 
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that f^l^n = 0, where ^ denotes the normal derivative with respect to the Riemann metric 



di^ \ dD _ 



(a^^) = fa,-i) ^, we have 



h{t, Xt) - h{0, Xo) = + 



where 



M,^ = h{t, Xt) - h{0, Xo) - ^ + ^) 



)ds 



is a martingale under P^, and 



A' 



t / g 

— + L ) h{s,X,)ds. 



(5.1) 



(5.2) 



On the other hand, applying Ito's formula 
h{t,Xt)-h{0,Xo) = 



d 1 A 



ds 2 ^ dx^dx^ 



h{s, Xs)ds 



it thus follows that 



and 



.7=1 



dx^ 



h{s,Xs)dMi 



(5.3) 



(5.4) 



Consider a solution u{x,t) to the initial boundary problem to the non-linear parabolic 
equation 

^ ^ d \ 

- L u + f{t,u,Vu) = 0, 



dt 

u{0, x) = ip{x), 
du{t 



0, 



t > 



(5.5) 



dD 



Then, by and 



hiT,XT) = h{t,Xt) + 

d „x 



9_ 

5s 



+ L h{s,Xs)ds 



^h{s,X,)dMi 
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together with the PDE fl5.5p we deduce that 

hit,Xt)-hiT,XT) = ^ fiT-s,h{s,Xs),Vhis,X,))ds 



[ ^Ks.Xs)dMl. (5.6) 



Let Yt = u{T — t, Xt) and = -£^h{t, Xt). Then the previous equation may be written 

as 

ZidMi (5.7) 



Yt-YT= [ f{T-s,Y,,Zs)ds-J2 f 
Jt Jt 



and Yt = u{0, Xt) = ip{Xt)- That is to say that Yt = u{T — t, Xt) solves the scalar BSDE 

d 

dYt = -fT{t,Yt,Zt)dt + J2 ZidMi, Yt = ^{Xt) (5.8) 

i=i 

where /t = /(T - t,y,z). 

For any fixed T > 0, let = {y^^ : t G [0, T] } be the unique solution to the BSDE (M 
on (f2, ^'^j ^t", P^). Since the solution to BSDE is unique, Y^ = u(T — t, Xt). In particular 
u(T, Xq) = Y^, and therefore 

I u{T,x)fi{dx) = E>'{Y^). (5.9) 

The above argument leading to the probabilistic representation fl5.9p can not be justified 
in the case that a = (a*-') is only Borel measurable or the boundary dD is only continuous, 
as in this case, {Xt)t>o is no longer a semimartingale, both (15.21) and (15. 4p no longer make 
sense. While, in this case, boundary problem (strong or weak solutions) to the non-linear 
PDE (15. 5p also need to be interpreted. On the other hand, the BSDE ( 15. Sp . which relies on 
only the martingale representation, still make sense, thus the representation theorems stated 
in §3 can be made as the definition of a solution to (15. 5p . This is the approach we will carry 
out. 

Consider the initial value problem of the following non-linear parabolic equation in a 
bounded domain D with a continuous boundary dD 



\dt 2 ^ dx^ ^ dx' 1 



u + f{t,u,Vu) = (5.10) 



5J = 

subject to the initial and boundary conditions 

d 

U{X,0) = Lp{x), ^^(^r 
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= for f > 

dD 



where a = (a^^) is Borel measurable, satisfying the uniform elhpticity condition: 

d d d 

A 5^ Iff < < A-^ Yl l^'l' ^(^^) ^ 

i=l i,j 2=1 

for some constant A > 0. 

Definition 5.1 The functional on ^i{D) defined by fi ^ {y(t,/i)o}, denoted by u(t, fj.), 
is called the stochastic solution of the initial and boundary problem of 1^5. 10\) . where for each 
t > and n G ^i{D), Y{t,ii) = {Ys)s<t is the unique solution to the BSDE 

d 

dYs = -fit - s, Ys, Zs)ds + ZidMi, Yt = ^{Xt) (5.11) 

on {n,^'^,^i',Xt,et,p>^). 

As a consequence we have 

Theorem 5.2 If (p is bounded and Borel measurable on D, and f is Lipschitz continuous, 
then there is a unique stochastic solution to the non-linear parabolic equation ( 15. 

We will study the regularity theory of the stochastic solutions in a separate paper. On the 
other hand we would like to derive an alternative probability representation of the stochastic 
solution. 

Let us apply the approach outlined in j2Tj. Let Yg = Ns — Vs where 



.-1 "''0 



Then V = (ys)se[o,t] is the unique solution to the functional differential equation 

Vs = r fit - r, Y{V)r, Z{V)r)dr, Vo = (5.12) 



where N{V)s = -E^'i^iXt) + Vtl^^, 

Y{V)s = E^^iipiXt) + Vt\^,^} - Vs 

for s G [0,t], and Z{V) is given as the density process of N{V) in the martingale represen- 
tation. In particular 

YiV)o = E^^{ipiXt) + Vt\^>;}. 
We therefore have the following 
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Theorem 5.3 Let ip be bounded and measurable. For t > 0, let V{t) be the unique solution 
to the functional differential equation liS.lS^) . Then the stochastic solution to the Neumann 
boundary problem of the non-linear PDF Ii5.10\) is given by 

u{t, ^) = {ip{Xt) + V{t)t} V// G yi(D). (5.13) 
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